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1. INTRODUCTION 



Energy spectra of complex quantum systems often exhibit a crossover from the 
Poisson statistics of uncorrected levels to the Wigner-Dyson statistics, with its 
characteristic level repulsion and rigidity. There is therefore some interest in 
designing and studying random matrix ensembles which could interpolate 
between these two extremes [1]. One such generalisation goes back to the work of 
Rosenzweig and Porter [2] on complex atomic spectra. They have studied, 
numerically, random matrices with enhanced (as compared with the standard 
ensembles) weight for the diagonal elements. A slight generalisation of the 
Rosenzweig-Porter matrix model leads to the following ensemble for N x N 
hermitian matrices : 

H = A + \y (1.1) 

N 

where A is a random diagonal matrix, i.e. 

A = diag(aiL a^^) (1.2) 

with some probability distribution p(a) for its statistically independent elements. 
V is a real or complex hermitian matrix, whose matrix elements are supposed to 
be independent, with a gaussian distribution of unit variance. There is no loss of 
generality, by considering A to be diagonal as long as we are interested only by 
correlations functions of the eigenvalues of H. 
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Different behaviour can be expected by varying the exponent a. If a = 2 / one 

expects GOE or GUE statistics whereas if a > 1 the statistics should be of the 

Poisson type. It seems that the only way to obtain new statistics is to choose a = 1. 

For technical reasons, up to now only the case of a complex hermitian V matrix 

has been treated. In an important work, which has inspired us, Brezin and 
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Hikami [3], abbreviated in the following BH, showed that when a = 2 / and for any 

fixed matrix A, the statistics is GUE. Pandey [4] has discussed the case a = 1 in 
great generality, using a Fokker-Planck equation for the probability distribution of 
eigenvalues. In the complex case, he could successfully use the fact that this 
equation is integrable and gave an expression for the two-point correlation 
function and the spectral form factor in a number of cases, including the one 
discussed here. 



Unfortunately, no derivation of the result was given. More recently, Guhr and 
Guhr and Miiller-Groeling [5] have derived a double-integral representation of 
the two-point correlation function, using a variant of the supersymmetric 
technique. In our opinion, however, this work rests on unproven prescription for 
obtaining the correlation functions from averages of advanced Green's functions 
only. 

Finally, Altland et al [6] using perturbation theory in A, obtained a closed 
expression for the two-point correlation function. However, since they kept the 
imaginary part of the energy non zero, before going to the large N limit their 
expression gives only an "overview" of the picture : it is limited to large distances 
and completely misses the fine structure of the correlation function at shorter 
distances. 

In the present paper, we look again at this problem (in the complex case), using 
BH strategy. We show first that the density of states v(E) is given by p, and that if 
we unfold the two-point correlation function, in the usual way in energy, and 

also in the coupling constant X replacing it by then we get a universal 

function, which coincides with that given by Pandey. Our approach gives more 
naturally the spectral form factor. Both quantities can be expressed as integrals of 
Bessel functions with imaginary arguments. Quite generally, we find that the 
model gives repulsion of levels at short distances and attraction of levels at long 
distances. The small coupling limit (close to Poisson) and the strong coupling 
limit (close to GUE) are described explicitly. In the small coupling limit, we 
recover the results of Leyvraz and Seligman [7] obtained by a perturbation 
technique. 



2. PROBABILITY DISTRIBUTION OF EIGENVALUES AND DENSITY OF STATES 



We consider an Hamiltonian H of the form 



H = A + V 



(2.1) 



A being a fixed matrix and V a matrix belonging to the unitary ensemble of 
probability distribution 



P (V) ~ exp 



^ tr 
2o 



(2.2) 



a being a parameter whose dependence on N will be fixed later. 



The probability distribution of H will be therefore 



P (H) ~ exp 



1 2 1 

— tr H + - tr AH 
2a o 



(2.3) 



In order to obtain the probability distribution W(?i^...Xj^) for the eigenvalues of 

H, one integrates over the unitary matrices U diagonalising H, by using the 
Harish-Chandra, Itzykson-Zuber identity [1] 



JdU 



exp 



- tr AmU^ 



A(X) A(a) 



det exp 



(2.4) 



the variables a^ being the eigenvalues of the matrices A and A(a) being the van 
der Monde determinant 



(2.5) 



c is a constant depending on N. 



If we use the identity 



. dX^ A(^) exp 



1 N 

i=i 



N 



= A (a) (2710)2 



(2.6) 



whose proof is given in the appendix, for completeness, we get for the normalised 
probability distribution W of the eigenvalues 



V 1 N7 



N 



r 1 A 



\2%o J 



1 A(X) 
A(a) 



exp 



det exp 



- a. X 



(2.7) 



In what follows, we will always consider average values of symmetric functions 
of the eigenvalues F(A,^, 'X^), with respect to the distribution W. In such a case 

the average value can be simply written as 



(F>w = 



f 1 A 



N 

2 1 



A(a) 



|d?ij...d?ij^ F(?i^...?Ij^) A(?i) exp 



(2.8) 



because the determinant of the matrix exp 



- a. X. 



can be simply replaced by 



N! exp 



Thus the full average of a symmetric function F over the ensemble of H matrices 
defined by Eq. (2.1) will be 



N 

F = jda,L da^, Hpl^i) (F 

i=l 



since we consider a model where all the a^ are statistically independent and 
identically distributed. 

Following BH, it appears more convenient to compute first the Fourier 
transforms of correlation functions, i.e. 



Ci(t) = E exp (it X^) 



C2(t, t') = £ exp (it + it' 



(2.9) 



Since from (2.6) follows that 



|dA,j L dA,j^ exp 



A(A,) exp 



(2.10) 



N 



= A (a + ba) (2 7t a) 2 exp 



we get 



Ci)w= exp-^ £ exp(itaj n 



1 + 



ita 



a,, - a. 



k "j 



and 



(2.11) 



kl 



1 + 



ito 



n 



1 + 



it'o 



X F (a^ , a,) 



(2.12) 



/ N a, - a, a, - a, + it'o- ita 

where F (aj, , a-^j = ^ ^ ^ 



a-^ - a^. - ita a-^ - aj^ + it' a 



(2.13) 



It remains to average over the a's. For this purpose we follow BH and use the 
following integral representations for (Cj)^ and (Cj)^. 



i/w 



exp 



f 2^ 

at 



2 , 



^ I ^ e'^p (it^) n 



ita J? 27ui 



i=k 



1+ ^ 



(2.14) 



and 



C2>w = -e'^P 



'^R 'r 



N 

j=i 



1 + 



ita 



z-a- 



1 + 



it' a 



z-a. 



F (z, z') 



(2.15) 
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being the contour formed by a rectangle of vertical size 2 e and horizontal one 

2R encircling the origin. One needs to choose R > sup a- , and in (2.15) one needs 

3 ' 



to impose that |t| > 2e and |t'| > 2e, in order to avoid the poles in the function 
F(z,z') defined by (2.13) 



F (z, z') 



z'-z - ia(t-t') z'- 



z'-z +it'o z'-z - ita 



(2.16) 



In this way, one gets for the Fourier transform of the density of states 



q(t) = exp-|t^ ^ f^exp(itz) 



f I ^ 

1 + iot < > 

^ z-a J 



N 



(2.17) 



where from now, < • > will designate an average with respect to the variable a, i.e. 



<^— > = [da p(a) — 
z-a 



z-a 



(2.18) 



Up to now, the validity of Eq. (2.17) is guaranteed only for a probability 
distributions of the a's of bounded support. However, since, as it is easily checked, 
we can take the limit R ^ <» in this formula, the contour becoming the 

contour F made of the two lines x ± ie. In this way the formula becomes valid 
for distribution p(a) of unbounded support. 

Equation (2.17) holds for any N. It is instructive to consider the N ^ oo limit, in 
order to exhibit the importance of the parameter a. 

We can rewrite Eq. (2.17) as 



^ = ^ ''"p 



ot 



r dz itz 
r^"^ j=o 



N 

I i J 



(iat) 



z-a 



(2.19) 



The j = term vanishes. The j = 1 term does not depend on N. The behaviour of 
higher terms depends on the value of the parameter Na. Since we are interested 
in the case where Na tends to zero, we see that all terms beyond j > 2 disappear so 
that 



lim — Ci (t) 



dz e 



itz 



27ui 



1 ita 

< > = <e > 

z-a 



(2.20) 



Hence the average density of states v(A,) is simply given by p (A,) [6] in our case. 



THE SPECTRAL FORM FACTOR 



We want now to compute the N = <» limit of the correlation function C(t, t') 
defined by 



NC (t, f) = C2 (t, f)- q (t) q (f) 



(3.1) 



From Eq. (2.14) and (2.15), we see that we can represent it in the form 



C (t, t') = [K, (t, t') + (t, t')] expf- I (t'+ t'') 



(3.2) 



with 



^1 = ^ii ^ K "'^ " p^^"/'^ p^^"'/"'^ 

r r 



and 



(3.3) 



■K, (t, t') = , ^ , f — I— exp[i(tz + t'z')l (z, z') 



1 1 

+ 



Z-Z'+T z'-Z + T' 



(3.4) 



where we used the notation 



T = ita 



t' = it'a 



g (z, z') = 1 + X a(z) + x' a(z') + ^fa(z')-a(z)] 

z — z 



(3.5) 
(3.6) 



with a (z) = < > 

z-a 



(3.7) 



and 



p(z/x) = 1 + X a(z) 
p(z'/x') = 1 + X' a(z') 



(3.8) 



From now on, we will take o = — ^ 

N 



(3.9) 



In order to get a finite limit for C, it is appropriate to choose the scale, defined by 



t = - + Ns 
2 

T 

t' = - - Ns 



(3.10) 



and keep T and s fixed when N tends to oo. We then choose e = — in the 
^ N 

contour F and represent the variables z and z' by 



z = X + 



Z = X 



2N 



y 



i5 
N 

16 



(3.11) 



2N N 



q and q' being equal to ± 1, depending on the branch of F on which the variables 
stay. We note that the function a(z) defined by Eq. (3.7) is such that 



a 



iSq 



with 



= Yq (^) + 



(x) = P jP^^ + i 7. q p(x) 
^ x-a 



(3.12) 



Since x = and x' = , at the order — ^, it follows that 

N N ' N 



p(z/T) = 1+ Y rx)+ 



p(zVt') = 



N '1 



1 - ^ Y,. (x) + 



N 1 



InJ 



(3.13) 



and if q = q' 



g(z, z') = 1 + O 



^ 1 ^ 



(3.14) 



but if q = -q' 



g(z, z') = 1- p(x) 



1 + 



sX'^iq 
qy - 2i5 



+ O 



^ 1 ^ 



(3.15) 



From these results, one can derive the desired asymptotic behaviour 



(t, t') = £ jdx p(x) e 



iTx -IkX sqp(x)^ 



where 



=dy ib(y-i5) 

e 

'In 



exp 



' ia ^ 



(3.16) 



(3.17) 



with F-, 



F. = 



47tp 



y + i 



Xsq _ g 



n2 



(3.18) 



and 



-4 2 
a = 7CA. s p 

b = 2sq 



(3.19) 



It remains to evaluate these integrals. One first notice that they vanish when b is 
negative. This implies that in Eq. (3.16), one can drop the summation on q and 
replace sq by |s| . 

L-^ can be evaluated in closed form and is given by 



Li = L (2Vab) 



(3.20) 



Ij (z) being the modified Bessel function given by the series 



2j+l 



L (z) = V 



(3.21) 



In evaluating L2 for b positive, one notices that when b is positive 

2 2, , 

A* s X* s 

— ^ - 5 = — — - 5 is positive because the formula for K- is valid under this 
2 2^ 2 

28 26 

constraint only, for N very large because of the condition |t| > — |t'| > — 
that we needed to impose in Eq. (2.15). One can then express L- in the form : 



L2 = 



j Ii(2Vab(t + l)) exp 



47up 



2 2 ^ 



(3.22) 



Grouping all these results we can finally write the correlation function C in the 
form 



C(t, t') = -|dxp(x) e'^'' S 



u = 



P(x) 



/A = Xp(x) 



(3.23) 



where the function S (u / A) is given by 



S(u/A) = - - Ij(y) exp 



00 , . 

+ ^ Y jdt (t^-l) Ii(Yt) exp 



,2 ,22 
-27uA u - A u 



2.2 2 ^ .2 
-t A u - 27iA u 



(3.24) 
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with Y = i8%A u 



(3.25) 



These expressions are our main result. We prove now that S is the unfolded 
spectral form factor. 



Indeed, by definition, C(t, t') is related to the eigenvalues correlation functions. 



Pi w = E^l'^-^) 



(3.26) 



by the equation 



C(t, t') = jdx jdx' e^^*^ ' '""'^ 1 [P2(x,x') - p^(x) p^(x')] 



(3.27) 



Let us make in this integral the change of variables (v(E) being the density of 
states) 



x = E + 



x' = E 



2Nv(E) 



2Nv(E) 



(3.28) 



and go from the variables (t, t') to (T, s) as defined by equation (3.9). 



Then we can see that we can express C as 



C(t,t') = -jdEv(E) e'™ jdr (E,r) exp 



s ^ 

ir 



V 



v(E) 



(3.29) 



where 



E + 



2Nv(E) 



2Nv(E) 



-P2 



r r 

E + 3TT^, E 



V 



2Nv(E) 2Nv(E) 



(3.30) 



is the unfolded (in energy) connected two-point correlation function. 



Notice that in our case v(E) = p(E). 



Our result for the asymptotic behaviour for C can therefore be interpreted as 
saying that if we unfold (E, r) in the coupling constant too, i.e. replace A, by 
A\v(E), then (E, r) tends to a function Y(r/A) independent of v. 

and 

jdr Y(r/A) exp (iru) = S(u/A) (3.31) 

is the spectral form factor. In fact, we found S(u) = S(|u|). We will now analyse in 
more detail the unfolded form factor S(u / A) and the connected correlation 
function Y(r / A) . 



4. TWO-POINT CORRELATION FUNCTION. STRONG AND WEAK COUPLING 
LIMITS 

Let us first give some properties of S and Y valid for any coupling A . 

Developing in series I^(Yt) in Eq. (3.24), and integrating over i, one can rewrite S 
in the form 



S = 



27tA^u 



k=2 



2n 



kiy) 



exp 



.2 2 ^ ,2 

-A u - 27uA u 



(4.1) 



where Ij^ is the modified Bessel function of order k. 



This expression shows that S is positive, a property possessed by the GUE form- 
factor. Since the correlation function Y is given by 



1 poo 

Y(r) = — cos ur S(u) du 

7t -"O 



(4.2) 



we see immediately that Y(o) > and Y"(o) < , which show that we have 
level repulsion at short distances. 

Consider now the long distance behaviour. From (4.1) one can see that for small u 



S(u) ~27iAu 



(4.3) 



This implies that |dr Y(r) = 



(4.4) 



to be contrasted with the sum rule 



jdr Y(r) = 1 



(4.5) 



in the GUE case, characteristic of the long range coulomb repulsion of the levels. 



Moreover, an integration by parts in (4.2) gives 



Y(r) = - 



2A 



- [cos ur S"(u) du 
7ur ° 



(4.6) 



This shows that we have an attraction of levels at long distances, and that 

1 

correlations decay very slowly, like . 

r 

In order to study the asymptotic behaviour of S and Y, it is helpful to look for 
another integral representation of these quantities. 

It can be easily checked from the series expansion of the Bessel function Ij^, that 
we have 



f— 1 
V u ) 



ik(y) = I 



Z>1 



dz 1 



A u ^ .2 
+ 27tA uz 



(4.7) 



Inserting this representation in the equation (4.1) for S, we can write it in the form 



S = I 



dz 1 1 



Z>1 



27ui z z-1 



u 



1 - 

^ 27uz" ) 



exp 



- aV(i-z - 27lA^u(l-z) 



(4.8) 



when — >1, we choose for the integration path on z a circle of radius ■\— ■ We 

27t V 27U 

can do the same choice when — <1, if we take into account the contribution of 

27U 

the pole at z = 1. In this way, one gets 



S(u/A)-S(u/oo)=-- jdy 



u 



1-y 



- + 2y Ji^ + 1 
27t ^ V27t 



exp 



-27uA u 



— + 2y J — + 1 

27U ^ V27t 



(4.9) 



where 



S(u/oo) = e 



I 27rj 



^ u ^ 
1- — 

I 27i; 



(4.10) 



is the GUE form factor. 



The corresponding equation for the correlation function Y is 



4 7 1 2yvx+l 1-y 

Y(r/A) - Y(r/oo) = -- [dx cos27trx \= 

^ ^ ^ ^ % { \ x + 2yVx+l 



exp 



-{2nAf X [x + 2Vx y + l]] 



and 



(4.11) 



Y(r/oo) 



^sinTir^ 
{ nr J 



(4.12) 



This expression coincides with that of Pandey [4] if we identify our A with half 
his. 
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This representation of Y is useful to analyse the large A limit, because A appears 
as the parameter of a Laplace transform. The dominant contributions to the 
integral will come therefore from the neighbourhood of x = and of x = 1 and 
y = -l. 

A careful analysis shows that the contribution from x = gives : 



5Y^ = -2 



A 



(2.A^)^ 



+ 



+ r 



(-1 



and the contribution from x = 1, y = - 1 gives 



CO +00 



= +2 

A % 



Re exp(27cir) jdy jdx exp 



r 2 \ 
y 

2 2 



X +y 



exp 



2 2\ . 471U r 

X +y + 1^ 

^ I V2 A. 



+ A 



-5/2 



Evaluating the integral and adding these contributions we find that 



Y(r) = 



2A 



+ 



+ r 



2{KTy 



cos 27tr exp 



-1 



+ A 



-5/2 



(4.13) 



It is important to note that this asymptotic expansion holds uniformly in r, i.e. the 

f -5/2 A -5/2 

symbol A means a term bounded by a constant times A , for all r. 



If we look at the large A limit, from the point of view of the form factor, we can 
see that the neighbourhood of the Heisenberg time is critical, in the sense that the 
following scaling limit exist. 



lim 27tA fS(u/A)-S(u/oo)] u= 27i + — 
A^oo Ay 



1 7 dt / 2\ 



(4.14) 



Let us comment briefly about the comparison of our result, for the large A limit 
(Eq. (4.11)) with previous work. We find agreement with the result of Guhr and 
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Miiller-Groeling [5] only in the regime r<<A . Our result is consistent with that 
of perturbation theory [6] if we assume that r»l and average over few level 
spacing. In this case we get 



Vp(r) = - - 

r + 



2A 



+ 



2 ' 2 2 

In r 



which corresponds to the first two terms, which are dominant in Eq. (14) of [6] (the 

2 

X in this reference should be identified with 2A ). 

It remains finally to examine the small A limit. For this purpose, it appears useful 
to have another integral representation of Y. 



d^S 

Using equations (4.7) and (4.8), one can see that — j can be expressed in terms of 

du 

Bessel functions 1^. Using recursion formula for these functions and the 
expression (4.6) for Y, one can express the correlation function in the form 



Y(r) = -— Tdu si 
r 



sinur 



exp 



2 2 

- A u -27tAu 



+ ■ 



2A 7 cosur 



r J u 



exp 



2 2 

- A u - 27uAu 



(4.15) 



This expression shows that if one keeps — fixed and let A tends to zero, one 
obtains a well defined quantity 



Y(r) = J du sinu — exp-u +0 (VAj 



A 



(4.16) 



which coincides with the perturbative result of Leyvraz and Seligman [7]. 
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Appendix 



We give here a short proof of the basic identity 2, used in the text. 



Let 



G(a) = (27io) 2 I dA,jL dX^ A(A, + a) exp 




where 



A(X + a) = n (^i 




N>i>i>l 



It is easily seen that G(a) is a totally anti symmetric polynomial in the a's of 
degree N-1 in each variable a^ . Therefore 

G(a) = C„ A(a) 

But since G(a) = A(a) when o = 0, we conclude that G(a) = A(a), for all a. 
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